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ABSTRACT 

We present a novel statistical analysis aimed at deriving the intrinsic shapes and 
magnetic field orientations of molecular clouds using dust emission and polarization 
observations by the Hertz polarimeter. Our observables are the aspect ratio of the 
projected plane-of-the-sky cloud image, and the angle between the mean direction of 
the plane-of-the-sky component of the magnetic field and the short axis of the cloud 
image. To overcome projection effects due to the unknown orientation of the line-of- 
sight, we combine observations from 24 clouds, assuming that line-of-sight orientations 
are random and all are equally probable. Through a weighted least-squares analysis, 
we find that the best-fit intrinsic cloud shape describing our sample is an oblate disk 
with only small degrees of triaxiality. The best-fit intrinsic magnetic field orientation is 
close to the direction of the shortest cloud axis, with small 24°) deviations toward 
the long/middle cloud axes. However, due to the small number of observed clouds, the 
power of our analysis to reject alternative configurations is limited. 



1 INTRODUCTION 



Far-infrared and sub-millimeter emission from molecular 
clouds appears polarized, presumably as a result of the align- 
ment of elongated dust grains with the cloud magnetic field 
(see, e.g., Curran & Chrysostomou 2007; Dotson et al. 2000, 
2009; Hildebrand et al. 2000; Hoang & Lazarian 2008; Lazar- 
ian 2003, 2007). Measurable degrees of polarization, at the 
few percent level, are typical in many interstellar clouds 
and cloud cores (Dotson et al. 2000, 2009). Additionally, 
theoretical advancements in understanding grain alignment 
(e.g., Bethell et al. 2007) indicate that even in relatively 
dense clouds, dust polarization traces the underlying mag- 
netic field structure. The ubiquity of magnetic fields and 
dust in the interstellar medium makes polarization observa- 
tions a powerful tool for distinguishing between theories of 
molecular cloud formation, support, and evolution. 

From a theoretical point of view, the orientation of the 
mean magnetic field in molecular clouds is closely tied to the 
dynamical importance of magnetic forces compared to grav- 
ity, random motions (turbulence), and thermal pressure. If 
magnetic fields are dynamically important, and are respon- 
sible for a significant fraction of the support of molecular 
clouds against gravity, then the mean magnetic field is pref- 
erentially oriented parallel to the shortest axis of the molec- 
ular cloud (e.g., Mouschovias 1978). This is a result of an 
increased support against gravity in the direction perpendic- 
ular to the magnetic field compared to the direction parallel 



to it. The cloud contracts more in the direction parallel to 
the magnetic field than in the direction perpendicular to the 
field. An additional result of magnetic support is that the 
intrinsic shapes of molecular clouds in this case resemble 
mostly oblate (although not necessarily axisymmetric) , flat- 
tened ellipsoids (one axis appreciably smaller than the other 
two). 

If magnetic fields are dynamically unimportant com- 
pared to turbulence, then turbulent motions dominate 
the internal dynamics of clouds, and the magnetic fields 
are dragged around by turbulent eddies (e.g., Ballesteros- 
Paredes, Vazquez-Semadeni & Scalo 1999). In this case, the 
mean magnetic field has a random orientation with respect 
to the molecular cloud principal axes. The shape distribu- 
tion for overdensities formed in a turbulent field is also ran- 
dom (Gammie et al. 2003). Molecular clouds forming out of 
magnetohydrodynamic turbulence in the weak-field regime 
therefore would be expected to have random shapes and 
random magnetic field orientations. 

A third possibility is that magnetic fields have a heli- 
cal configuration and thread prolate (filamentary) molecular 
clouds (Fiege & Pudritz 2000a). The most common out- 
come of such configurations is polarization patterns that 
may contain 90° flips of the polarization vector (Fiege & 
Pudritz 2000b). Additional ideas for molecular cloud shapes 
come from non-magnetic calculations; for example, flnite, 
self-gravitating gaseous sheets have been shown to collapse 
to filamentary structures with mass concentrations close to 
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the edges or ends of the filaments (e.g., Burkert & Hartmann 
2004; Hartmann & Burkert 2007). 

If it were possible to determine the mean orientation 
of the ordered magnetic field in molecular clouds, then im- 
portant constraints could be placed on theories of molecu- 
lar cloud formation and support. However, such a task is 
not straightforward, due to various projection effects which 
prohibit us from knowing either the true orientation of the 
magnetic field or the principal axes of the cloud ellipsoid. 
Polarization measurements can only determine the orienta- 
tion of the magnetic field on the plane-of-the-sky (POS), 
and, similarly, only the shape of the POS projection of the 
cloud can be measured. As a result, the magnetic field may 
be oriented, for example, along the shortest cloud axis, but 
its projection on the sky may appear closest to the major 
axis of the projected cloud ellipse. These difficulties were 
explicitly demonstrated by Basu (2000). 

An effective way to overcome these difficulties and use 
polarization measurements to constrain theoretical models 
for the orientation of the mean magnetic field in molecular 
clouds is through a statistical treatment. Assuming that the 
orientation of the clouds themselves with respect to our LOS 
is random, we can assess the likelihood of different under- 
lying distributions of molecular cloud shapes, and of orien- 
tations of the mean magnetic field in molecular clouds. In 
this work, we present a new method for such a treatment. 
We study a sample of 24 molecular clouds with measured 
apparent elongations and field orientations. By applying a 
weighted least-squares analysis, we derive the most probable 
intrinsic shapes and most probable magnetic field orienta- 
tions for the parent population that the 24 clouds sample. 

This paper is organized as follows. The observations and 
analysis used to derive elongations and apparent field orien- 
tations in the sample clouds are described in !j2l The for- 
malism used in our analysis, including a short discussion of 
projection effects and a description of our statistical analy- 
sis, are discussed in ^ (a detailed mathematical treatment 
of projection effects and the statistical analysis are given in 
Appendices 1X1 andlBlrespectivelv). Our results are presented 
in ^ and discussed in ^ 



2 OBSERVATIONS 

The Hertz polarimeter (Schleuning et al. 1997; Dowell et al. 
1998) deployed at the Caltech Submillimeter Observatory 
has been used to map the flux density and polarization of 
molecular clouds on scales of several arcminutes at a wave- 
length of 350 /im. In a complete sample of 56 Galactic clouds, 
32 contain 2 or more polarization measurements satisfying 
the criterion P > 3ap or better, where P is the measured po- 
larization amplitude and ap is its measurement uncertainty^ 
The clouds are listed in Table [T] object positions and com- 
plete maps of flux density and polarization can be found in 
Dotson et al. (2009). 
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Figure 1. Mean polarization and cloud shape in DR21. Polariza- 
tion measurements (small black line-segments) are plotted with 
their lengths proportional to P (scale at lower-left); only P > 3(Tp 
measurements are shown. The mean polarization direction {E- 
vector) and its dispersion are shown as the red line and the arcs 
at its end. The cloud principal axes are drawn as blue lines; their 
relative lengths are given by the aspect ratio q = 0.26. The lines 
corresponding to the mean polarization and cloud widths inter- 
sect at the cloud center {xg^yo) and are drawn with arbitrary 
absolute lengths. The gray arc labeled "A" denotes the difference 
between the mean polarization angle and the cloud position angle. 
Coordinate offsets from 20'^39™in, 42°19'31" (J2000). Contours 
are 10, 20, 90% of the peak flux density of 820 Jy/beam (at 
(Ao, AS) = (0', 3f3)). The shaded circle at lower-right is the 20" 
FWHM Hertz beam. 



2.1 Polarization 

Within each cloud, we limit the data points used to cal- 
culate mean polarizations to those satisfying the condition 
P > 3(Tp. The best estimate of the POS-projected mag- 



^ In counting 32 clouds, we have counted separately two compo- 
nents of OMC-3 with distinct cores separated by more than 2', 
with no Hertz measurement of flux in between. 
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Figure 2. Same as Figure [T] but for the cloud W49A. As this 
cloud is nearly circular (aspect ratio q = 0.87) we plot only the 
long cloud principal axis for clarity. Coordinate offsets are from 
ig^'lO^lSfe, 9°6'17" (J2000). Contours are 10, 20, 90% of the 
peak flux density of 730 Jy/beam. 



netic field direction is that of the measured polarization vec- 
tor (rotated by 90°). The amplitude of the polarization P 
yields some information about the inclination of the field to 
the line-of-sight (LOS). However, P is also dependent on a 
number of other factors which are independent of the field, 
including the dust grain shape and alignment efficiency, and 
the magnetic field structure along the LOS and within the 
telescope beam (e.g., Draine & Lee 1985). Therefore, our 
analysis of the magnetic field will involve only the polariza- 
tion position angle. 

The mean polarization position angle (x) is defined by 
averaging Stokes parameters of unit magnitude such that 



(X) 
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(cos2x> = ;^ 5Zcos2xi, 
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and the Xi's are the measured polarization angles at each 
position i in the cloud. The sums are over all A'^ points in 
the map where the measured polarization P > 3ap. 

We calculate the standard deviation of the mean (x) 
under the assumption of zero intrinsic dispersion (i.e., under 
the assumption that all of the dispersion in the observed x* 
is due to observational errors) as 

1/2 



1 

iV 



(4) 



where the ai are the uncertainties on the measurements of 
Xi- The measurement uncertainty is limited to Hertz's sys- 
tematic uncertainties such that a(^-^^ > 2° (Dowell et al. 



1998). In our case, the dispersion in the observed values of 
Xi in individual clouds is substantial, so we also calcinate 
the dispersion about the mean, a^, given by the standard 
deviation of the measurements 

1/2 



iV 



(5) 



In using this equation we have accounted for the 180-degree 
periodicity in the polarization angles. 

Each cloud's mean polarization amplitude is given by 



(P) 



(6) 



where I, Q, and U are the measured Stokes parameters. 
The mean polarizations, angles, and uncertainties are all re- 
ported in Table [1] The {P)'s reported in Table [1] have been 
de-biased using the statistical measurement uncertainty on 
the mean (e.g., Vaillancourt 2006). However, the measure- 
ment uncertainties in the table are not those used to de-bias, 
but are limited to Hertz's systematic uncertainties such that 
> 0.2% (DoweU et al. 1998). 

2.2 Cloud Shapes 

In the case of a symmetric ellipsoidal cloud, the first- and 
second-moments of the flux density describe the cloud's cen- 
ter and width. The location of the cloud's center {xo,yo) is 
given by the first-moments: 

J2i F{xi,yi) 



yo = 



T,iF{xi,y^) 



(7) 
(8) 



where the sum is over all points i in the flux density map. 
No cuts are applied to the flux density data beyond those 
already made in the Dotson et al. (2009) archive. (The P > 
3ap criterion is not applied.) 

The second moments compose the 3 elements of a sym- 
metric 2x2 matrix 



^xx ^xy 
Ixy lyy 

with elements 

Ixx = '^{y^-yQf F{x,,yi), 



'Y{x^ ~ xo){y^ - yo) F[xi,yi), 



hy = 'Yix^- xof F{x^, 



(9) 

(10) 

(11) 
(12) 



The principal moments of the cloud are simply the eigen- 
values, Ji and I2, of I. Defining Ji > I2, the aspect ratio 
is given by g = ^/h/Ii and the position angle of the long 
cloud axis is determined from the eigenvector corresponding 
to the eigenvalue J2. 

The first and second moments are well-defined and ex- 
press the distribution of the flux density in the cloud even 
when the shape is irregular and cannot be well approximated 
by an ellipse. Therefore, we use these moments to define the 
aspect ratio and position angle of the principal axes for all 
clouds in our sample, with the caveat that we have ignored 
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any higher-order moments which more precisely characterize 
the cloud shape. We note that some bias can be introduced 
through the finite size of each map. 

Examples of the mean polarization and cloud shape cal- 
culations are shown in Figure[T]for the most elongated cloud, 
and Figure[2]for a nearly circular cloud. The measured cloud 
parameters are given in Table [1] The apparent cloud aspect 
ratio q and the absolute value of the angle between mean 
field and short apparent axis, A, are plotted in Fig. [3] Data 
shown as open circles in Fig. [3] and denoted by h in Ta- 
ble [T] are not included in further analysis due to their large 
(> 30°) dispersion in A. 

Error bars on A correspond to the quadrature sum of 
cr(^), the systematic angle uncertainty (2 degrees; Dowell 
et al. 1998), and the systematic polarization uncertainty 
(Axsys). We estimate the last quantity using the relation 



90 



90° 0.2% 

TT (P) ' 



(13) 



where 0.2% is Hertz's systematic polarization uncertainty. 
We note that equation (|13l) assumes measurements of P and 
X follow normal distributions (equivalent to the assumption 
P ^ CTp). While this is clearly not true for all measurements 
in Table [TJ this fact has minimal effect on the apparent size 
of the error bars and has no effect on our subsequent analysis 
(as the uncertainties are not used)0 

Furthermore, these angle uncertainty estimates do not 
include any uncertainties in the estimate of the cloud shape, 
but include only the uncertainties on the measurement of the 
polarization angle. A visual comparison of Figures [T] and [2] 
should make it apparent that the uncertainty in the cloud 
orientation is smallest when q = and must increase as the 
cloud becomes more circular {q — 1). An estimate of this 
uncertainty for all clouds in our study is beyond the scope 
of this paper. Neglecting these uncertainties has no effect 
on our conclusions as this information is not used in our 
subsequent data analysis. 

Data shown as open circles in Fig. [3] and denoted by h 
in Table [1] are not included in further analysis due to their 
large (> 30°) dispersion in A. 



3 ANALYSIS 

3.1 Projection Effects 

Let us consider a triaxial ellipsoid model molecular cloud 
with semi-axes a> b > c, and axial ratios 



(middle-to-long axis ratio) and 



(14) 



(15) 



(short-to-long axis ratio). When observed, the image of the 
molecular cloud appears on the POS as an ellipse. The ob- 
servable quantity related to the cloud shape is the ellipse 
aspect ratio, q. Let us additionally assume that the molecu- 
lar cloud is threaded by a magnetic field, the mean direction 

^ Further study of polarization angle uncertainties at low P/o-p 
can be found in Naghizadeh-Khouei & Clarke (1993) 
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Figure 3. Elongations and angles between the magnetic field 
and the minor cloud axis for the 32 clouds in our sample. Note 
that the elongation q is defined so that a very elongated cloud 
corresponds to g — > and a circular cloud to g = 1. For angles, 
A = corresponds to a projected magnetic field aligned with 
the minor axis of the projected cloud. Each error bar represents 
the quadrature sum of the statistical uncertainty of the cloud's 
mean angle uncertainty (ct^.,^^), the systematic angle uncertainty, 
and the systematic polarization uncertainty (see text for precise 
definitions of these terms). Data shown as open circles are not 
included in further analysis due to their large dispersion in A (all 
data are given in Table [T}. Dotted horizontal lines are drawn at 
the boundaries of the allowed range 0° < A < 90°. The diagonal 
A/90° = q is shown as a solid line. 
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Figure 4. Solid black lines: native cloud ellipsoid coordinate sys- 
tem {x,y,z). Dashed blue lines: observation coordinate system 
{x' , y' , z'). The LOS is along the z' axis, and the POS is the x'-y' 
plane. The y' axis represents the direction of the projection of 
the shortest ellipsoid axis onto the POS. The directions of the 
magnetic field and of its projection onto the x-y plane are shown 
in red. 



of which forms a polar angle 9b with the short axis of the 
cloud ellipsoid, and an azimuthal angle (j>B with the long 
axis of the cloud ellipsoid (see Fig. |4| . Through polarimetry 
observations, only the direction of the projection of the field 
on the POS can be measured, and the associated observable 
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Table 1. Mean Polarization and Cloud Shape Parameters 



Mean Polarization Parameters {E-vector)^ Cloud Shape 



Object 


Object,^ 


No. 




o-(p> 


(x> 






Aspect Ratio 


Angle_£ 




No. 


Name 


vectors 


(%) 


(%) 


(degrees) 


(degrees) 


(degrees) 




(degrees) 


(degrees) 


1 


W3'^ 


51 


0.6 


0.2 


65 


2 


32 


0.83 


—80 


35 


2 


NGC 1333 


18 


1.6 


0.2 


88 


2 


19 


0.68 


—45 


48 


3 


IRAS 05327-0457 


18 


3.9 


0.4 


—40 


2 


17 


0.92 


—37 


3 


4 


OMC-1 


437 


2.0 


0.2 


25 


2 


26 


0.51 


5 


20 


5 


OMC-2'' 


26 


0.7 


0.2 


—39 


4 


33 


0.44 


14 


53 


6 


OMC-3 MMS 6_£ 


43 


2.1 


0.2 


-42 


2 


12 


0.74 


-51 


9 


7 


OMC-3 MMS 9-'' 


21 


1.2 


0.2 


-62 


2 


20 


0.95 


11 


73 


8 


TV T A 

OMC-4 


16 


1.7 


0.2 


— 14 


2 


13 


0.84 


57 


71 


9 


NGC 2024 


54 


0.7 


0.2 


—23 


2 


24 


0.65 


— 15 


8 


10 


NCG 2068 LBS17 


3 


8.3 


1.5 


81 


5 


16 


0.98 


—50 


49 


11 


NCG 2068 LBSIO 


43 


3.8 


0.2 


—51 


2 


20 


0.93 


38 


89 


12 


NGC 2071 


6 


0.6 


0.2 


—27 


4 


25 


0.97 


— 13 


14 


13 


Mon R2" 


49 


0.7 


0.2 


31 


2 


32 


0.82 


47 


16 


14 


GGD12 


17 


1.2 


0.2 


89 


2 


16 


0.94 


—60 


31 


15 


S269 


7 


2.9 


0.3 


30 


3 


13 


0.82 


— 18 


48 


16 


AFGL 961 


5 


1.9 


0.7 


— 10 


4 


44 


0.86 


—37 


27 


17 


Mon OBI 12_£ 


26 


2.0 


0.2 


—37 


2 


17 


0.76 


27 


64 


18 


NGC 2264 


18 


0.6 


0.2 


77 


2 


16 


0.84 


—47 


56 


19 


p Oph 


41 


1.6 


0.2 


-21 


2 


13 


0.76 


7 


29 


20 


IRAS 16293 


7 


0.6 


0.2 


90 


3 


25 


0.67 


-65 


26 


21 


NGC 6334 


8 


0.3 


0.2 


90 


3 


34 


0.93 


-53 


37 


22 


NGC 6334 a"" 


49 


1.3 


0.2 


69 


2 


18 


0.89 


- 3 


71 


23 


NGC 63341 


54 


0.9 


0.2 


41 


2 


19 


0.51 


11 


30 


24 


W33 C*" 


29 


0.4 


0.2 


48 


2 


37 


0.88 


82 


35 


25 


W33 A 


19 


0.9 


0.2 


-53 


2 


13 


0.90 


-49 


4 


26 


M17 


127 


0.9 


0.2 


-13 


2 


27 


0.59 


-10 


3 


27 


W43-MM1 


4 


1.5 


0.3 


9 


4 


22 


0.88 


34 


25 


28 


G34'' 


44 


0.5 


0.2 


-66 


2 


41 


0.94 


72 


42 


29 


W49 A 


32 


0.6 


0.2 


59 


2 


29 


0.87 


-80 


41 


30 


W51 A 


109 


0.5 


0.2 


44 


2 


28 


0.71 


-78 


58 


31 


W75 N'' 


9 


0.3 


0.2 


44 


3 


47 


0.98 


76 


33 


32 


DR21 


142 


1.2 


0.2 


5 


2 


20 


0.26 
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" Mean polarization parameters defined in equations JTJ - (|6j. Means are calculated using only the P > 3o"p vectors shown in "No. 
vectors" column. Angle x is measured east of north. 

Object names follow the convention in Dotson et al. (2009); list is ordered with approximately increasing Right Ascension. 

Dispersion about mean polarization angle given by equation JSj. 

Direction of long cloud axis, measured east of north 
^ Angle between mean polarization angle and long cloud axis. This is equivalent to the angle between the inferred magnetic field direction 
and the short cloud axis. 
/ Part of OMC-3 in Dotson et al. (2009) 
9 a.k.a. IRAS 06382-f0939 

^ Objects with angle dispersion larger than 30 degrees are not included in the angle analysis. 



quantity is the angle A between the projected field and the 
minor axis of the POS cloud ellipse. 

The observables q and A can be calculated as a function 
of ^, Ob, 4>b and of the orientation angles of the observer's 
LOS, 6 and (j). This calculation is described in detail in Ap- 
pendix [X] 

However, the orientation of the LOS is unknown — in 
the absence of biases, the LOS orientation can be treated 
as uniformly distributed among different directions with re- 
spect to the observer. As a result, it is not generally possible 
to de-project the intrinsic cloud shape and magnetic field 
orientation for any single object. Instead, a statistical treat- 



ment must be used: the preferred intrinsic cloud shapes and 
magnetic field orientations in nature can only be obtained 
by observing a large number of clouds, and comparing the 
distribution of observables with expectations based on dif- 
ferent intrinsic shapes and field orientations. In this work, 
we employ such a statistical analysis described below. 



3.2 Statistical Analysis 

As a first proof-of-concept, we test for the presence of possi- 
ble degeneracies in the distribution of the observables (g. A) 
which may limit the value of this analysis. The simplest 
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tests we can perform to verify that the data of Fig. |3] have 
nontrivial information content derive from the immediate 
observation that almost all datapoints lie below the diag- 
onal extending from (0.0,0) to (1.0,90). We evaluate the 
probability to obtain this configuration by chance. 

First, we consider the case in which both aspect ratios 
and angles are randomly drawn from uniform distributions. 
Under this assumption, the desired probabilities can be ob- 
tained analytically: the diagonal splits the parameter space 
in two parts of equal area. Therefore, the probability for any 
single point to lie above or below the diagonal is 0.5, and the 
outcome of multiple draws obeys the binomial distribution 
with probability for a positive outcome equal to 0.5. If we 
consider all 32 datapoints (including the ones that did not 
survive our quality cuts), then we have 2/32 points above 
the diagonal. The probability to obtain the observed out- 
come or one which is even more biased towards points below 
the diagonal is the sum of the probabilities to obtain 0/32, 
1/32 or 2/32 points above the diagonal: 

P(# above diagonal < 2) = ^ ^.^ 0.5'0.5^^"' = 1.2x10"^ 

(16) 

If we consider only the 24 datapoints that survived our qual- 
ity cuts, then we have 1/24 points above the diagonal, and 
the associated probability to have obtained such a result by 
chance is 

P(# above diagonal < 1) = X] f ^•'^) 0.5*0.5^""* = 1.2x10"'^ 

(17) 

with the result being numerically approximately equal as in 
the previous case. This configuration is extremely unlikely to 
have been obtained by chance from a uniformly distributed 
parameter space. 

Next, we consider the case in which the aspect ratio dis- 
tribution in nature is identical to the observed aspect ratio 
distribution, but angles are randomly drawn from a uniform 
distribution. For each observed aspect ratio, we randomly 
select an angle A, and we calculate the probability that the 
number of points below the diagonal is < 2 for the 32 dat- 
apoints and < 1 for the 24 datapoints surviving the cuts. 
We find that this probability is equal to 0.98% and 1.3% 
respectively. Although more likely than before, this scenario 
also has a very small probability of occurring. 

A more sophisticated version of these qualitative ar- 
guments can be made starting from characteristic limiting 
cases for the distribution of intrinsic cloud features (shapes 
and magnetic field orientations) . For each such case we com- 
pute, by performing "mock observations" along randomly 
selected and uniformly distributed LOS, the joint proba- 
bility density function (PDF) of our two observables, the 
elongation q and the projected field - short ellipse axis angle 
A. The results are shown in Fig. [SI overplotted with the 24 
observed points that survive the A dispersion cut described 
in !j2l The parameters of the intrinsic shape and field ori- 
entation distributions that lead to each PDF are given in 
Table (2] These parameters are: in the case of the intrinsic 
shape distribtion, the cloud ellipsoid axial ratios (Cmax, ^max) 
of the most probable shape (the values that maximize the 
PDF of axial ratios); and in the case of the field orientation 
distribution, the most probable cosine of the field orientation 



polar angle cos 0s, max, and the standard deviation of the co- 
sine of the field polar angle, (j(cos0_g). The exact functional 
forms of the shape and field orientation distributions used 
are discussed in Appendices IB II and IB2I 

The qualitative distribution of the data points on the 
g-A plane forms a triangular shape, with the lower-right part 
of the plot (almost circular cloud images, small angles of the 
magnetic field with the short ellipse axis) more populated 
than the upper left part of the plot (very elongated cloud im- 
ages, large angles of the magnetic field with the short ellipse 
axis). In general, there are many fewer very elongated POS 
cloud ellipses than circular-looking ellipses, which qualita- 
tively indicates oblate intrinsic shapes (since there are many 
more LOS that will yield a circularly-looking POS ellipse for 
a disk-like cloud than for a cigar-like cloud). Additionally, 
large angles with the short ellipse axis are only encountered 
for almost circular cloud images (giving rise to the triangu- 
lar distribution of the data on the observables plane) , which 
qualitatively hints to an orientation of the magnetic field 
at small angles with the short cloud axis (so that when the 
cloud is viewed edge-on, the observations yield an elongated 
POS ellipse and a small POS field angle with the short el- 
lipse axis, cf. Fig. [1] while when viewed face-on the cloud 
looks circular and a larger variety of angles of the POS field 
with the short ellipse axis are possible, cf. Fig. [2]). 

The different PDFs plotted in Fig. [5] demonstrate how 
the behavior of the PDF of the observables responds to 
changes in the underlying distributions of shapes and field 
orientation. For oblate shapes, the contours of the joint PDF 
for A and q for difi'erent field orientations are shown in panels 
a-d. As expected, when the field is oriented parallel to the 
short cloud axis, the X-q PDF acquires a roughly triangu- 
lar shape preferentially populating the high-g low-A corner 
(panel a), while for a field oriented perpendicular to the 
short cloud axis the PDF is also roughly triangular, how- 
ever now preferentially populating the high-q high-A corner 
(panel b). For the intermediate situation of a field oriented 
at 45° from the short cloud axis, the PDF forms a "spike" 
at the mid-g mid- A part of the plot (panel c). Finally, for 
a random distribution of field orientations (panel d), the 
contours are parallel to the A axis, with low g values pre- 
ferred. PDFs resulting from prolate intrinsic cloud shapes 
are shown in panels e and f, with panel e corresponding to 
a magnetic field oriented along the long axis (filamentary 
clouds formed along the magnetic field), and panel e corre- 
sponding to a magnetic field oriented along the short axis 
(filamentary clouds perpendicular to the magnetic field) . In 
both cases, low values of g are preferred (a situation not 
seen in the data). The PDF of panel e preferenially popu- 
lates the high-A part of the plane, while the PDF of panel f 
preferentially populates the low-A region. 

We can therefore see that there is enough qualitative 
and quantitative variation between these PDFs so that, with 
a sufficiently large dataset, the effect of the random orien- 
tation angles can be overcome and we can draw conclusions 
regarding the underlying distributions of cloud shapes and 
magnetic field orientations in nature. 

In order to determine which such intrinsic distributions 
of cloud parameters fit our data best, we bin the A-g param- 
eter space in 5 X 5 bins. For each bin, we count the number 
of observed data points that fall within the bin limits. We 
parameterize the intrinsic shapes distribution as described 
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Figure 5. Contour plots of the joint probability density function of the observables q and A obtained by convolving different intrinsic 
distributions of cloud shapes and magnetic field orientations with random LOS. The solid, dashed, and dotted lines correspond to the 
1,2, and 3(t contours, respectively. The blue data points correspond to our 24 data points surviving the quality cuts. Panel a: oblate 
clouds, magnetic field preferentially oriented parallel to short cloud axis; panel b: oblate clouds, magnetic field preferentially oriented 
parallel to the long cloud axis; panel c: oblate clouds, magnetic field preferentially oriented at a 45 degree angle with the short cloud axis; 
panel d: oblate clouds, uniformly distributed magnetic field orientations; panel e: prolate clouds, magnetic field preferentially oriented 
parallel to long cloud axis; paned f: prolate clouds, magnetic field preferentially oriented parallel to short cloud axis. The parameters of 
the intrinisc distributions resulting in each PDF are given in Table [2] 



Fig. [5] 
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oblate, B || short axis 
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0.999 
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0.19 


oblate, B ± short axis 


c 


0.999 
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0.04 


oblate, B at 45° with short axis 


d 


1.00 


0.62 


uniform 
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oblate, B random 


e 


0.29 


0.08 


0.19 


0.19 


prolate, B± short axis 


f 


0.29 


0.08 


0.91 


0.16 


prolate, B|| short axis 



Table 2. Parameters of the intrinsic distributions giving rise to the PDF in each panel of Fig.[5l 



in Appendix IB II using a bi-parametric joint distribution in 
C, and ^ of non-zero spread, uniquely defined by the val- 
ues of C, and ^ where the probability density becomes max- 
imum, C^max and ^max. We parameterize the intrinsic distri- 
bution of magnetic fields orientations as described in Ap- 
pendix [B2l using a bi-parametric distributiorl3 in cos 6s, 

^ Since a likelihood analysis on the shapes alone (see discussion 



uniquely defined by the value of cos 6b where the proba- 



in Appendix IBll l indicates that the distribution of cloud shapes 
is stronly peaked on the C = 1 axis, implying that the long and 
middle axes of the cloud ellipsoid are very close to being equal, 
the <j)B angles are physically degenerate, and the only quantity of 
interest determining the orientation of the magnetic field is the 
angle with the shortest axis of the cloud ellipsoid, Ob- 
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bility density becomes maximum, cosSsmax, and the stan- 
dard deviation of the distribution, a(cosS_B). Note that we 
wish to know the distribution of cos 6b rather than that of 
6b, since the quantity of interest is the fraction of LOS per 
sohd angle in any given direction 9b, which is expressed by 
cos 6b (see also discussion in Appendix IB)) . Assuming uni- 
formly distributed LOS, we determine the joint PDF for q, 
A, and the expected counts for 24 observations in each bin, 
as well as the spread in the expected counts. We perform a 
weighted least-squares analysis to find the best-fit parame- 
ters Cmax, Cmax,cos 6'_Bmax, cr(cos 6b), usiug the iuverse square 
of the spread in the expected counts as our weight. The de- 
tails of this statistical analysis are discussed in Appendix 
[B] Because of the nonlinearity of the problem, there are 
no applicable analytic solutions to the minimization of the 
weighted sum of squared residuals. For this reason, we have 
used a well-tested, Monte-Carlo-based approach, the "Sim- 
ulated Annealing" algorithm (Corana et al. 1987) to sample 
the parameter space and identifying the location of the ab- 
solute minimum of the weighted sum of squared residuals. 



4 RESULTS 

The results of our weighted least-squares analysis are shown 
in Fig. |6] The upper panel shows contours of weighted 
squared residuals, ^weighted, as defined in Appendix |B] on 
the Cmax-^max plauc, with COS 6'smax and a (cos 6lfl) fixed 
at their best-fit values. The lower panel shows contours of 
^weighted on the cosSB,max-o-(cos6'_g) plane, with i^max and 
^max fixed at their best-fit values. The allowable values of 
a(cos&s) vary between (5-function of B-field orientations) 
and 0.289 (uniform distribution of B-field orientations with 
all orientations equally probable, see Appendix|B} . The color 
scale corresponds to values of Swoightod from 9.7 to 25, with 
the contours spaced by factors of 1.1 (so the suppressed z- 
axis in this plot which is visualized by the contours is in log- 
arithmic scale). The location of the minimum is indicated 
in each case by the white point within the back region. The 
outermost contour corresponds to the Swoightod value that is 
typically yielded by 24 observations drawn from the best-fit 
intrinsic distributions due to random fluctuations, as cali- 
brated by Monte-Carlo simulations. 

As expected from the qualitative arguments outlined 
above, oblate cloud shapes are preferred. In the case of the 
magnetic fleld orientation, there is a well-defined broad min- 
imum in the weighted least-squares analysis, corresponding 
to a clear preference for a small angle from the short cloud 
ellipsoid axis and a moderate spread of angle about that. 
However, given the small number of data, the power of the 
test to reject models is limited, and a large fraction of the 
cos 6 Bma.x-o' {cos 6b) plauc (excludiug, however, distributions 
strongly peaked at large angles from the shortest cloud axis) 
could be consistent with the presently available data. 

The best-fit values of the parameters are shown in the 
first line of Table (2] and the resulting PDF of our observ- 
ables, q and A, is shown in panel a of Fig.[Sl The best- fit dis- 
tribution of intrinsic shapes (axial ratios of the model cloud 
ellipsoids) and magnetic field orientation angles are shown 
in the upper and lower panels of Fig. [T] respectively. The 
most probable shape is a moderately thick oblate disk. The 
best-fit shapes distribution is very strongly peaked close to 
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Figure 6. Contours of the weighted sum of squared residuals, 
•S'woightod, for the (Cmax,€max) parameters (upper panel) and the 
{cos 8b max: ""(cos parameters (lower panel). The upper axis 
in the lower panel shows the most-probable orientation in de- 
grees from the short cloud ellipsoid axis, corresponding to the 
value of cos(9s_niax shown in the lower axis. The two parameters 
not shown in each plot are kept fixed at their least-squares val- 
ues. The color scale corresponds to values of Swoightod from 9.7 
to 25, with the contours spaced by factors of 1.1 (so the sup- 
pressed z-axis in this plot which is visualized by the contours is 
in logarithmic scale). The minimum is indicated in each case by 
the white point within the back region. The final contour corre- 
sponds to the 5wcightod value that it typically yielded by 24 ob- 
servations drawn from the best-fit intrinsic distributions due to 
random fiuctuations, as calibrated by Monte-Carlo simulations. 
The solid line in the upper panel separates the domain of allowed 
values ((^0 ^ ?o) from the rest of the plane, while the dotted 
line separates mostly oblate [c^max > 0.5(1 + ^max)] from mostly 
prolate [C^max < 0.5(1 + ^max)] most-probable ellipsoids. 



« 1.0 (so clouds are likely to have very small degrees of tri- 
axiality). On the other hand, the distribution is very spread 
out in ^, implying that many different disk thicknesses are 
possible. Although the most probable shape is a relatively 
thick disk, this does not mean that all clouds are thick disks. 
Much thinner disks [closer to ((^, ^) = (1,0)], as well as very 
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Figure 7. Best-fit intrinsic distributions of cloud axial ratios (up- 
per panel) and of magnetic field orientation angles (lower panel). 
In the upper panel, the red line extending from (CiO=(lil) to 
((^, ^)=(0.5,0) separates the p(CiO = plane in two sections: el- 
lipsoids with axis ratios to the left of this line are mostly prolate, 
and ellipsoids to the right of the line are mostly oblate. Infinitesi- 
mally thin disks have ((^, 5)=(1,0); infinitesimally thin cigars have 
(f,^)=(0,0); and perfect spheres have ((^, 5)=(1,1). 

thick clouds [closer to (i^, ^)=(1,1), labeled "sphere" on the 
p(C,,£,) ~ plane], are also frequently encountered although 
not equally common in this distribution. 

The most probable magnetic field orientation has a 
small offset {9bo ~ 24°) from the short ellipsoid axis, and 
the standard deviation in cos Ob of ~ 0.21. These results 
are robust against the removal of any outliers from our 
dataset of 24 observations. For example, when repeating 
the analysis without our most elongated cloud (DR 21), the 
best-fit parameters of the shape distribution change from 

(Cmax,Cmax) = (0.999,0.58) tO (Cmax.fmax) = (0.986,0.64), 



while the best-fit parameters of the B — field orien- 
tations distribution change from [cos Svmax, o-(cos Ss)] — 
(0.91,0.21) to [cos6'vmax,a(cos6'B)] = (0.86,0.19) - a change 
smaller than 10% in any one of the parameters, while the 
qualitative behavior of the underlying distributions remains 
unchanged. This result is not surprising: in the upper-left 
panel of Fig. [2]DR 21 is outside the 2a contour of the best- 
fit joint PDF for q and A, and barely within the 3a contour. 

As indicated by the broadness of the minimum of 
^weighted as a fuuctiou of our fitted parameters, the best-fit 
intrinsic distributions of shapes and field orientations are not 
unique in their ability to yield q-X PDFs that are acceptable 
representations of the data used in this analysis. However, 
both the qualitative arguments presented in i)3.2l as well as 
the location of the minimum of Swcightcd indicate that oblate 
shapes and mean magnetic field orientations with small de- 
viations from the shortest axis of the cloud are preferred. 



5 DISCUSSION 

In this paper, we have presented a new method for assess- 
ing quantitatively the intrinsic shapes and the orientations 
of the mean, ordered component of the magnetic field of 
molecular clouds, using observations of the apparent aspect 
ratios q and angles A between the apparent mean magnetic 
field and the apparent short axis of the POS cloud images. 
Under the assumption that the LOS towards the observed 
clouds have random orientations with respect to the princi- 
pal cloud axes, we have explored different statistical meth- 
ods to evaluate the consistency of various models for the 
intrinsic shapes and intrinsic magnetic field orientations of 
clouds with our data. We have used data from 24 molecular 
clouds obtained through 350 ^im observations with the Hertz 
polarimeter (Dotson et al. 2009). 

Based on our data sample we can exclude certain sim- 
ple scenarios with high confidence. A scenario in which both 
cloud aspect ratios and magnetic field orientation angles are 
randomly drawn from uniform distributions can be excluded 
at the 10~^ level. A less restrictive scenario, in which the 
distribution of cloud aspect ratios is identical to the one ob- 
served in the data but magnetic field orientations are drawn 
randomly from a uniform distribution, is also excluded at 
the 10"^ level. 

In order to examine more general cases, we have em- 
ployed a weighted least-squares analysis to derive the best- 
fit intrinsic shapes and intrinsic magnetic field orientations 
of the clouds in our sample. We have found that the most 
probable intrinsic shape is a thick (short to long axis ra- 
tio ~ 0.6) oblate disk with a negligible degree of triaxiality 
(middle to long axis ratio ~ 0.99); the best-fit distribution 
of cloud thicknesses was found to be broad, so both thin 
and thick clouds are frequently encountered under this dis- 
tribution. The most likely orientation of the magnetic field 
is close to the shortest axis of the cloud (~ 24° offset toward 
the middle/long axis). The best-fit distribution of magnetic 
field orientations is clearly peaked around this value, but 
also features long tails. We have found that the best-fit dis- 
tribution of shapes and orientations, when convolved with 
a random LOS distribution yields a distribution of observ- 
ables which is in good agreement with the data. Our results 
are robust against the removal of any single point in our 
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dataset (for example, if we repeat the analysis without our 
most elongated cloud, DR 21, the change in our best-fit pa- 
rameters is smaller than 10%). 

These results are in agreement with the qualitative 
trend scon in the data, where clouds with small apparent 
aspect ratio (apparently elongated clouds) have small an- 
gles between the mean projected magnetic field direction 
and the short axis of the cloud ellipse; clouds which are ap- 
parently circular on the other hand can have a large range 
of angles between mean magnetic field and short axis. Fi- 
nally, apparently elongated clouds with large angles between 
mean field direction and short cloud axis are not seen in the 
data. As far as the distribution of apparent cloud shapes is 
concerned, more clouds are seen to have large aspect ratios 
(and are thus close to circular) while only a few clouds are 
apparently very elongated. 

Indeed, clouds which are intrinsically disks with the 
magnetic field direction close to that of the shortest cloud 
axis would exhibit similar properties when viewed through 
random LOS. More LOS yield a disk seen almost face-on 
rather than edge-on, so most clouds would look almost cir- 
cular and only a few would appear significantly elongated 
(in contrast to prolate clouds, most of which would appear 
elongated, and randomly shaped clouds, which would show 
no preference in apparent aspect ratio). In addition, when 
clouds are seen almost edge-on and appear elongated, the 
magnetic field will be aligned with the short axis of the pro- 
jection of the cloud on the POS. In contrast, when clouds arc 
seen face-on, the magnetic field projection could from any 
angle with respect to the short and long axes of the cloud 
ellipse, which in reality are, in this case, projections of the 
cloud ellipsoid middle and long axes. 

Despite the preference of the data for such a scenario, 
and the clear potential of this test for discriminating be- 
tween scenarios for intrinsic magnetic field orientations in 
molecular clouds and between theories for molecular cloud 
dynamics, wo have shown quantitatively that a larger num- 
ber of observations is needed for the test to conclusively 
reject alternative configurations. With the present, limited 
number of observed clouds, the distributions of LOS, mag- 
netic field orientations, and shapes, arc still sparsely sam- 
pled, and as a result random fluctuations in this sampling 
allows for a large number of parameters to yield acceptable 
representations of the observed q and A. However, our results 
explicitly demonstrate that there is no instrinsic degeneracy 
in the distributions of q and A yielded by different classes of 
cloud shapes and magnetic field orientations in nature, so 
a large number of q, A observations should allow conclusive 
tests of different underlying distributions. 

Our analysis has the advantage that it is not a priori 
tied to any theory or prediction regarding the dynamical 
processes in molecular clouds, but rather allows the data 
to pick freely the part of the shapes/orientations parameter 
space that best fits the observations. We have additionally 
tested that the location of the best-fit Cmax, Cmax, cos^smax, 
and the spread of the orientations distribution are robust 
with respect to changes in the assumed functional form or 
the details of the analysis (for example, performing a joint 
analysis for the shapes and orientations, or performing an 
analysis for the shapes alone). Finally, observational uncer- 
tainties in the measured quantities are expected to have a 
limited effect in our analysis, as the binning of the (g. A) 



parameter space involves bins that are typically of the same 
order as or wider than such uncertainties. 

The apparent orientation of the magnetic field in molec- 
ular clouds and cloud cores has been studied in the past by 
various authors. Kane et al. (1993), using 1.3 mm polariza- 
tion measurements, claimed a strong preference for align- 
ment of the polarized emission with with the structure in 
deconvolved IRAS maps. Glenn et al. (1999) used a sample 
of 7 elongated cloud cores selected to have a polarization 
detection greater than 3(t at either 800 /xm or 1.3 mm. They 
found that the orientation of the apparent magnetic field 
was random with respect to the apparent cloud axes. Vallee 
& Bastion (1999), using 760 /im observations, found the ap- 
parent magnetic field direction in molecular cloud cores (in- 
tensity peaks within clouds) to be parallel to the apparent 
minor axis in 3 out 10 cases. These studies all focused on po- 
larization vectors measured at intensity peaks, contrary to 
our own survey, whore polarization is extensively mapped 
throughout the clouds in our sample. The intrinsic shapes 
of molecular clouds were discussed by Kerton et al. 2003. 
Based on the Five College Radio Astronomy Observatory 
Outer Galaxy Survey in CO emission, they studied a sample 
of 15,000 clouds. They concluded that the intrinsic shapes 
of these objects are best described as intermediate between 
near-oblate and near-prolate ellipsoids. 

Our results have important implications concerning 
the dynamical importance of magnetic fields in molecular 
clouds. If magnetic fields arc dynamically unimportant com- 
pared to turbulent motions, the fields are expected to be 
carried around with turbulent eddies in the cloud, and show 
little correlation with the principal cloud axes. However, our 
data show clear indications that such a correlation exists. 
Our results are also unfavorable for the scenario of helical 
magnetic fields threading prolate clouds. Not only axe the 
predicted 90° flips in the polarization vectors not observed 
(Dotson et al. 2009), but clouds in our sample arc found to 
be oblate, rather than prolate. This result is also consistent 
with the argument against prolate cloud cores, which are 
also predicted in this scenario (Fiege & Pudritz 2000c), but 
are not observed in nature (Jones et al. 2001, Jones & Basu 
2002, Tassis 2007). 

If, on the other hand, magnetic fields are dynamically 
important, cloud shapes should be close to oblate disks, and 
the magnetic field should be closely aligned with the short 
axis of the cloud ellipsoid. This scenario seems to be the 
one which yields the best agreement with our data. Intrinsic 
cloud shapes are indeed very strongly peaked around oblate 
disks. The best-fit magnetic field orientation is close to the 
shortest cloud ellipsoid axis, with small offsets toward the 
middle/long axes. However, more data are required before a 
more robust quantitative statement can be made with regard 
to both the details of the best-fit underlying magnetic fields 
orientation distribution as well as the confidence with which 
alternative scenarios can be rejected. 

In addition to support from our statistical treatment, 
this picture is also consistent with observations of individ- 
ual well-studied systems where projection effects are mini- 
mal. In these cases, the hourglass morphology of the mag- 
netic field, characteristic of dynamically important magnetic 
fields, is seen in very different scales, from clouds to dense 
cores (e.g., Schleuning 1998, in the case of the cloud OMC- 
1; Girart, Rao & Marrone 2006, in the case of protostellar 
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core NGC1333; Lai et al. 2002, in the case of massive core 
NGC2024). 

We conclude by stressing that although our analysis has 
yielded best-fit distributions that seem to prefer this portion 
of the parameter space, the rejection power of our test, es- 
pecially in the case of magnetic field orientations, is limited 
given the present sample size. Stronger constraints on the 
intrinsic statistics of magnetic field orientations in nature 
will require polarimetry observations in a larger number of 
clouds. 



and ^ — c/a, so 1 > > 5. Then, the POS isophotes of an 
ellipsoid molecular cloud, the shape of which is described by 
Eq. HAip . will be coaxial ellipses of aspect ratio 



where 

A-: 



A + C- ^{A^cy + B^ 



A + c + ^iA-cy^ + B^ 



sin (h + 



+ 



(A2) 



(A3) 
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APPENDIX A: PROJECTION EFFECTS 

Let us consider a triaxial ellipsoid model molecular cloud 
of semi-axes a > b > c, and a native system of coordinates 
{x, y, z) centered on the cloud and aligned with its axes. The 
triaxial ellipsoid surface in this system obeys 



2 2 
a2 62 



+ 



(Al) 



The X, y, and z axes are thus oriented along the longest, 
middle, and shortest ellipsoid axes respectively. Let us addi- 
tionally define a second, observation system of coordinates, 
{x',y',z'). In this system, the LOS is along the 2' axis, and 
consequently the x'-y' plane is the POS. The x' axis is in 
the x-y plane, which implies that a;' _L 2. The y' axis thus 
represents the projection onto the POS of the 2 axis, and 
consequently of the shortest axis of the cloud ellipsoid. The 
orientation between the native and observation systems of 
coordinates is given by the angles 9 and <j) between the 2 
and 2' axes and the x and x' axes respectively (see Fig. |4ll . 

When observed, the image of the molecular cloud ap- 
pears on the POS as an ellipse. The properties of this pro- 
jected, observed ellipse can be calculated from the properties 
of the cloud ellipsoid and the orientation of the observer's 
LO^(Binney 1985). We define the axial ratios as ^ — b/a 



B = cos6lsin20 ( 1 - l^j 1^ . 



C = 



+ cos 



1 



(A4) 



(A5) 



The orientation of the axes of the POS cloud ellipse 
with respect to the POS axes (a::' and y') is given by the 
angle ip, defined as 



1 / B 

_arctan(-^ 



(A6) 



In Eq. HA6p . tp can be the angle between either the minor 
or the major axis of the POS ellipse with the y' axis (the 
sum of the two angles is always equal to 90°, and ip, which 
can take values between and 45°, is always the smaller of 
the two). Which one of the two angles is represented by ip 
in each case is determined by the sign of the quantity 



{A - C) cos 2V' + B sin 27/> . 



(A7) 



If M < 0, V represents the angle between the minor ellipse 
axis and y'. If u > 0, V represents the angle between the 
major axis and y' . The sign of ip indicates the quadrant in 
which the relevant axis (major or minor, depending on the 
sign of u) resides. If 1/' < the relevant axis is located in 
the first/third quadrant, while if t/) > the relevant axis is 
located in the second/fourth quadrant. 

Additionally, let us consider that the cloud is threaded 
by a large-scale ordered magnetic field B, which has some 
orientation with respect to the cloud principal axes charac- 
terized by a set of angles {8b,4'b)', here, 9b is the angle of 
the magnetic field with the shortest cloud axis 2 and 0s is 
the angle of its projection onto the plane {x,y) defined by 
the directions of the middle and longest axes (see Fig.|4]). We 
are interested in the orientation of the POS projection of B, 
since this is what can be measured through polarimetry ob- 
servations. The components of B in the {x,y,z) coordinate 
system are 



(sm Ob cos (ps 
sin 9b sin b 
cos 9 b 



(A8) 



Performing the rotation to the {x' , y' , z') coordinate system 
we find that the x' and y' components of B (the POS com- 
ponents) are 

B^i = B [~ sin (/)sin 9b cos (jis 4- cos </> sin 6^5 sin (Pb] (A9) 



* Note that the projection we are referring to here is not a "cut" 
of the cloud ellipsoid along some plane, but is rather the surface 



brightness along the observer's LOS (and hence it represents an 
integral of the cloud density along the LOS) 
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and 

Byi = B[-~ COS (f) cos 9 sin cos i^b 

— sin (ji cos 6 sin 6*5 sin + sin 6 cos ^b] 



(AlO) 



The angle u) between the POS projection of B and the POS 
projection of the shortest cloud axis (the y' axis) is given by 



tantj ■ 



sin sin 9b cos (^s — cos <j) sm 9b sin 



cos9sm9 B (cos(j)cos(l) b + sinijisini/js) — sinScosS b 

(All) 

Note that, contrary to tp, the angle ui is, by definition, posi- 
tive in the first-third quadrant, and negative in the second- 
fourth quadrant, where x' and y' have opposite signs. 

The angle that can be determined observationally, how- 
ever, is not uj, but rather the magnitude of the angle A be- 
tween the POS projection of B and one axis (for definiteness 
we will use the minor axis) of the POS cloud ellipse. The re- 
lation between A, ui, and is different, depending on the 
value of u (see Fig. IA1[I : 



A : 



90° 

190° 



|90° - IV'- 

■ IV" + ^1 1 , 



u<0 
u>0 



(A12) 



We have therefore shown that the two observable quan- 
tities q and A can be calculated using the formalism de- 
scribed above, provided that we know: 

(a) the intrinsic shape of the cloud ellipsoid (i.e., its axial 
ratios ( and ^); 

(b) the orientation of the ordered magnetic field with re- 
spect to the native coordinate system (i.e., the angles 9b 
and 0s); and 

(c) the orientation of the LOS with respect to the native 
coordinate system (i.e., the angles 9 and (p). 

Note that the angles 9b and (pB only enter the calculation 
of A, while q is not dependent on them. Because the LOS 
is random, A is not fixed, even if the true magnetic field di- 
rection is fixed with respect to the true cloud ellipsoid axes, 
and all molecular clouds have the exact same intrinsic shape 
(see e.g., Basu 2000). Rather, there is a distribution of such 
angles, which can be calculated analytically, given our hy- 
pothesis for the orientation of the true mean magnetic field 
in the clouds, as well as some knowledge of the distribution 
of intrinsic cloud shapes in nature. 



APPENDIX B: STATISTICAL ANALYSIS 

In this work, we use our set of observations of the aspect 
ratio q of the POS cloud ellipse and the magnitude of the 
angle A between the POS projection of the magnetic field 
and the minor axis of the POS cloud ellipse to constrain the 
intrinsic shapes and magnetic field orientations in molecular 
clouds. Unfortunately, it is not possible to derive directly 
the intrinsic shape and the orientation of the magnetic field 
in each cloud, as this would require knowledge of the LOS 
orientation in each case. Instead, we employ a statistical 
analysis, assuming that the distribution of LOS orientations 
is uniform with respect to all directions (same fraction of 
LOS per unit solid angle in all directions). In practice, this 
means that the distribution of cos 9 is uniform with values 
between —1 and 1, and the distribution of (p is also uniform 
with values between and 2tt. 
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Figure Al. POS angles between the POS projection of the mag- 
netic field (red line), the POS projection of the shortest cloud 
ellipsoid axis (the y' axis), and the projected cloud ellipse axes 
(blue lines). Upper panel: u < (ip represents the angle between 
y' and the minor ellipse axis). Lower panel: u > {ip represents 
the angle between y' and the major ellipse axis). The magnitude 
of the angle A is the observable angle between the POS cloud el- 
lipse minor axis and the POS projection of the magnetic field. In 
this case, ip and ui are both positive, and |i/> + < 90°. 



Under this assumption, we will perform a weighted 
least-squares analysis to deduce the best-fit probability den- 
sity function of our observables, q and A, and the associated 
distributions of intrinsic shapes and orientations of the mag- 
netic field with respect to the principal cloud ellipsoid axes 
that generate it. 

Because the magnetic field orientation does not enter 
the calculation of q, it is also possible to use only the ob- 
servations of q to derive the best-fit distribution of the in- 
trinsic molecular cloud shapes (an analysis similar to that 
of Tassis 2007 in the case of molecular cloud cores). The 
results obtained by this g— only analysis indicate that the 
distribution of intrinsic shapes is strongly peaked on the 
= 1 axis (largest and middle axes equal) so that all az- 
imuthal magnetic field angles (pB are physically equivalent. 
For this reason, we will treat the orientation angle tp>B as 
random and uniformly distributed, and we will only at- 
tempt to determine the best-fit distribution of cos 9b, de- 
termining the most likely deviation of the magnetic field 
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orientation from the shortest cloud axis. Note that the rel- 
evant quantity is cos 9b, rather than 6b, as we are not 
only interested in the most probable angle, but also on 
the spread of such angles: one of the theoretical cases we 
would like to test is that where all magnetic field orienta- 
tions have an equal probability of occurring. This translates 
to a uniform distribution of magnetic field orientations per 
solid angle around the cloud center, which in turn is math- 
ematically expressed as a uniform distribution in cos 6b- 
Thus, the interpretation of a best-fit distribution in cos 6b 
is intuitively straight-forward: the standard deviation cor- 
responding to a (5-function in preferred field orientation is 
o"(cos^fl) = 0, while the standard deviation correspond- 
ing to a uniform orientation distribution is (T(cos^_g) = 
V{(cose»s)2) - {cos6'i3)2 = ^1/3- 1/4 = 0.289. Moder- 
ately spread distributions will have a standard deviation in 
between. 

Our primary statistical analysis is performed under the 
following assumptions. 

i. Molecular clouds can be described as triaxial ellipsoids 
of the form of Eq. (|A1|I . 

ii. The axial ratios and ^ of molecular cloud ellipsoids 
have an intrinsic distribution in nature, of the form de- 
scribed in Appendix IBll with two free parameters ^max and 
^max (the C and ^ which maximize the distribution and cor- 
respond to the most probable shape). 

iii. <j)B has a uniform distribution between and 27r; 
cos 6b has an intrinsic distribution of the form described 
in Appendix IB2I with two free parameters cos Ss, max, and 
a {cos 6 b) (corresponding to the most probable orientation 
and the standard deviation of the distribution of cos 9b)- 

iv. The orientation angles of the LOS with respect to the 
native coordinate system, 6 and </!>, are random and all 
equally probable: the distribution of cos 6 is uniform with 
values between —1 and 1, and the distribution of (j) is uni- 
form with values between and 360°. 

Based on these assumptions, we calculate, for every 
set of distribution parameters ("max, Cmax, cos Ss, max, and 
a {cos 6 b), the expectation value of the number of observed 
points in each bin on the q-X plane, if 24 observations are 
performed with random LOS; we also calculate the spread of 
the observed number of points. We do so using the following 
Monte-Carlo procedure. 

(i) We randomly draw a pair of and a cos 6b from 
their respective distributions, and a 0s from a uniform dis- 
tribution between and 2n- 

(ii) We randomly draw a LOS orientation, by drawing a 
pair oi 6, cj> from a uniform probability distribution (equal 
probability for any value of cos S or in the intervals (0, 1) 
and (0, 27r), respectively). 

(iii) We use the values of ^, ^, 9, and (j> to calculate q, u, 
and %l) through Equations (fA2| - (|A7|) 

(iv) We use the values of 9b, 4>b, 6, and </!> to calculate cu 
through Eq. (jXTTJ. 

(v) We use the values of cj, ^, and u to calculate A through 
Eq. dHU. 

(vi) We repeat steps (i)-(v) 24 times (as many as our 
observed clouds). For every bin in the q-\ plane, we calculate 
the number of observations that fall within its limits. 

(vu) We repeat step (vi) 10,000 times. Using the 10,000 



mock sets of 24 observations, we calculate for each bin i 
the mean number of observations, (A^,;) (which has to be 
< (Ni) < 24), and the sample variance of Ni, cr^j = 
(Nf) — (Ni)'^- Because we only use 10,000 mock observations 
and we use the inverse of a%i as a weight in our least-squares 
analysis, for the cases wherein our 10,000 trials we obtain 
a%i = we use instead a floor value of a%i = 10"* 

We finally obtain the best-fit set of parameters ("max, 
Cmax, cos Ss, max, and (t(cos^_b) by minimizing the weighted 
squared residuals function, 

o f^^ o I- a \\ Adata,i — (Ni) 
Owcightcdl^Qmax, ^max, cos y_B.max, 0"(,COS y_BjJ — 2 

(Bl) 



Bl Modified Lognormal Distribution of Intrinsic 
Cloud Shapes 

We wish to use a distribution of finite width in each of 
the ^ axes which is, however, smooth and has a relatively 
sharp maximum. Such a distribution is a modified lognor- 
mal, which we construct in the following way: we first seek 
an appropriate minimal transformation which will transform 
the domain of ( and ^ from (0, 1) to {—oo, oo); we then take 
the transformed variable to follow a gaussian distribution. 

An appropriate transformation is ^ x and ^ —* y 
where 
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It is easy to see that 
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If then the joint probability distribution function (PDF 
of X and y is 



p{x,y) = 



■ exp 



(3 



- xo 



2(7? 



exp 



{y - 2/0 ) 



2ai 



(B4) 



then the joint PDF of ( and ^ will be 



p(C,e) = 



2-KaxO-y 



X exp 
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2al 
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This is in principle a tetra-parametric distribution with 
parameters a:;o, yo, fx, and Oy- However, to reduce the num- 
ber of parameters entering the problem and the associated 
computation time required for the analysis, while at the 
same time allowing the expected amount of variation in 
molecular cloud shapes, we fix the value of and Gy to 
0.2 and 0.7 respectively, determined from a likelihood anal- 
ysi^ of the aspect ratios q alone. 

^ For details on the methodology of this shapes-only likelihood 
analysis see Tassis 2007 
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The values C, — ("max and ^ = ^max that maximize ^) 
are given by the system of equations 



dp/d£, = 
dv/d(, = 

or, equivalently, 



(B6) 
(B7) 



= oli In -r^^ - xa ) (C - C) 



1-C 



yo) (1 - C)C + (1 - 2C + 0C<^''^' , 

(B8) 
(B9) 



With Ox and (j^ fixed, (^max and ^max can be used aher- 
natively as the free parameters of p(CjO- 



B2 Modified Lognormal Distribution of Intrinsic 
Magnetic Field Orientations 

The distribution of cosOb follows a modified lognormal sim- 
ilar to the one discussed in Appendix [Bl] let the variable w 
follow a Gaussian distribution with parameters wo,cjw, and 
let cos Ob be related to w through the transformation 



w = In ■ 



cos Ob 



1 — cos 9b 
In this way, the PDF of cos 6b is 



p (cos 6s) = 



27ro- 



exp 



1 — COS p 
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cos 6b (1 — cos 6 



(Bll) 

The value of cos 9b where p{cos9b) is maximized, which we 
call cos 6smax, is obtained by requiring that dp/ d cos 9b ~ 0, 
or, equivalently, by solving the equation 



cos 9b 2 , „ 2 a 

In 7— Wo — + 2cr„ cos 9b 

1 — cos 9b 



0. 



(B12) 



The standard deviation of this distribution is obtained 
through 



(j(cos 9b) 
with 



{cos9b — {cos 9b}) p{cos9B)dcos6B , 

cos 9^ —0 

(B13) 



(cosOb) = cos6Bp{cos9B)dcos6B ■ 

J COS R — 



(B14) 



This is a bi-parametric distribution, with free parame- 
ters Wo and a^, or, equivalently, cos^snmx and cr(cos6_g). 
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